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Abstract 

We demonstrate that supersymmetric decays, as typified by the pre- 
dictions of several GUT-scale boundary condition choices, do not prevent 
detection of H^A^,H+H-, at a 1 TeV - 4TeV e+e" or /i+zx" col- 

lider operating at anticipated luminosity. For much of parameter space the 
relative branching ratios for various SUSY and non-SUSY decays can be 
measured with sufficient accuracy that different GUT-scale boundary condi- 
tion choices can be distinguished from one another at a very high confidence 
level. 

1 Introduction 

The minimal supersymmetric model (MSSM) is widely regarded as the most 
attractive extension of the Standard Model (SM). The approximate unification of 
coupling constants that occurs in the MSSM at an energy scale of a few times 
10^^ GeV suggests the appropriateness of treating the MSSM in the context of 
a grand unified (GUT) model, in which the supersymmetry breaking parameters 
have simple universal values at the unification scale, M[/. The GUT framework 
is especially compelling in that electroweak symmetry breaking (EWSB) is easily 
induced at a scale ~ mz as the soft mass-squared of the Higgs field that couples to 
the top quark is driven to small (sometimes negative) values by the associated large 
Yukawa coupling during evolution to low energy scales. Thus, it is important to 
consider the implications of GUT scenarios for the detection of the Higgs bosons 
of the MSSM and to determine the extent to which (and strategies by which) 
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Higgs boson decay branching fractions can be measured with accuracy sufficient 
to constrain GUT models. 

The Higgs sector of the MSSM is reviewed in Refs. [Q, ^. The MSSM con- 
tains exactly two Higgs doublets, leading to two CP-even Higgs bosons {h^ and 
with rriho < rriHo), one CP-odd Higgs boson {A^) and a charged Higgs pair 
(H^). Crucial parameters for the Higgs sector are m^o and tan/? (the ratio of the 
vacuum expectation values for the neutral Higgs fields that give mass to up-type 
and down- type quarks, respectively). A fundamentally important GUT result is 
that essentially all models with proper EWSB require m^o > 200 GeV, with much 
larger values being common. This result has many important implications: 

• The will be very SM-like, and, at fixed tan/5, will have a mass near 
the upper bound predicted by including (two-loop/RGE-improved) radiative 
corrections as computed for the known value of rrit and the values for stop- 
squark masses and mixing predicted by the GUT. For all scenarios considered 
(even those with my^o well above a TeV), 171^0 is below ~ 130 GeV and, as 
reviewed in Ref. will be discovered with relative ease at both the LHC 
and any e^e~ or fj^~^fi~ collider with ^/s > 500 GeV. However, because the 

will be very SM-like, it will be quite difficult to establish on the basis 
of precision measurements that it is the MSSM and not the SM Higgs, 
especially if m^o > 300 - 400 GeV §. 

• The H^, and H"^ will be approximately degenerate in mass and will 
decouple from the vector boson sector. The coupling of the A^ to bb [tt\ is 
given by 75 times —gmb/{2mw)t8.nP [—gmt/ {2m w) cot (3]. For large m^o, 
the couplings of the if" asymptote to i times these same coefficients. The 

— > tb coupling is proportional to ig / {\/2mw){mi,PRtBxil3 + mtPi cot P). 

• In most GUT scenarios, the high masses predicted for the if ° and A^ imply 
that decays to pairs of supersymmetric particles will be important when 
tan f] is not large and ti decays are not kinematically allowed. For small 
to moderate tan/? and m//o,m^o > 2mj, tt is the dominant mode unless 
the mass of the lightest stop squark, ti, is small enough that decays to titi 
are kinematically allowed. (This does not happen in the GUT models we 
consider.) When tan/? is large, the enhanced bb coupling of the A'^ and 

imply that bb decays will become dominant, even when SUSY and/or 
tt decay modes are allowed. In the case of the if^, SUSY decays always 
compete with the larger tb decay mode since mH± > rrit + mj, for the GUT 
scenarios considered. (In the GUT models we consider, tibi decays are not 
kinematically allowed.) 

• For rriAo ^ 200 GeV it is entirely possible that none of these heavy Higgs 
bosons could be detected at the LHC (see the review of Ref. |^), even as- 
suming the absence of SUSY decays. In terms of the (m^o, tan (3) parameter 
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space plane, heavy Higgs discovery is not possible once m^o > 200 GeV if 
tan /? > 3 and if tan (3 lies below an upper limit that increases with increasing 
m^o (reaching tan /? ~ 15 by m^o ~ 500 GeV, for example). More than likely, 
the tan /5 < 3 discovery region would be diminished after including the SUSY 
decays of the and that are predicted to be important. Detection of 
the and via such SUSY decays at the LHC appears to be very difficult 
except in rather special situations. 

• The only large rate production modes for these heavy Higgs bosons at an 
e^e~ or ix^ iJi~ collider will be —>■ H^A^ and Z* — > H^H~ . These modes 
are kinematically limited to my^o ~ m^/o ~ wi//± < \/s/2. In particular, 
at a first e"'"e~ collider with = 500 GeV and L = 50 fb~ observation is 
restricted to roughly < 220 — 230 GeV, implying that detection would not be 
possible in most GUT scenarios 

• Although single H'^ and A^ production is significant at a 77 collider facility 
for masses < O.S^i, i.e. about 400 GeV at a ^/s = 500 GeV e~^e~ colhder, 
backgrounds are such that very high luminosities are required for discovery 
— L > 200 fb~^ is required when either SUSY decays are significant or 
tan /? is large. 

In combination, these results imply that if", A", and detection may very 
well require employing the Z* —>■ H^A^ and Z* —>■ H^H~ production modes at an 
e'^e~ or colhder with y/s substantially above 500 GeV. Even if the if" and 

A^ are observed at the LHC, studying their decays and couplings would be much 
simpler in the pair modes. Various aspects of Higgs pair production are discussed 
in Ref. 0, which appeared as we were completing the present work. 

Our first goal is to determine the luminosity required to guarantee observability 
of the Z* H^A^, H^H~ modes regardless of the SUSY-GUT decay scenario. 
We will consider collider energies of 1 TeV and 4 TeV (the latter being actively 
considered |^ for ^'^ ^~ colliders), with integrated luminosities up to 200 fb^^ and 
1000 fb~^, respectively. Our second goal will be to develop strategies for organizing 
the rates observed for physically distinct final states so as to yield information 
regarding the relative branching fractions of different types of decay modes, and to 
assess the extent to which such information can determine the GUT scenario and 
its parameters given the expected experimental errors. 

We find that if the integrated luminosity at 1 TeV (4 TeV) is close to 200 fb~^ 
(1000 fb"^) then detection of the i/^o and H+R- pair production processes will 
be possible over almost all of the kinematicaly allowed parameter space in the 
models we consider, but that significant reductions in these luminosities will imply 
gaps in parameter space coverage. A measurement of the mass m^o ~ mj^^o ~ 
iriHi already provides critical constraints on the GUT model. The correlation 
between this mass and the masses of the charginos, neutralinos, and/or sleptons 
(as measured in direct production) determines the GUT scale boundary conditions 
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(provided there is universality for the standard soft-SUSY-breaking parameters) 
and a fairly unique location in the parameter space of the GUT model so singled 
out. In particular, tan/3 is determined. Assuming full luminosity, the relative 
Higgs branching fractions can be used to cross check the consistency of the GUT 
model and confirm the parameter space location with substantial precision. For 
example, the relative branching fractions for the H^, A^, to decay to SUSY 
pair particle states vs. Standard Model pair states provide a surprisingly accurate 
determination of tan/? given a measured value for m^o. This tan/3 value must 
agree with that determined from the masses. Other relative branching fractions 
provide complementary information that can be used to further constrain the GUT 
model, and can provide a determination of the sign of the Higgs superfield mixing 
parameter. Thus, a relatively thorough study of the full Higgs sector of the MSSM 
will be possible and will provide consistency checks and constraints that could 
single out the correct GUT model. 

The organization of the paper is as follows. In the next section, we describe the 
six GUT models that we consider, and delineate the allowed parameter space for 
each. Contours of constant Higgs boson, neutralino and chargino masses arc given 
within the allowed parameter space, and Higgs boson decay branching fractions are 
illustrated. In Section 3, we demonstrate that, for expected integrated luminosities 
at e'^e~ or n'^iJ.~ colliders, detection of Higgs pair production will be possible in 
final state modes where both Higgs bosons decay to final states containing only b or 
t quarks, even though the branching fractions for such final states are decreased due 
to competition from the SUSY decay channels. Event rate contours as a function 
of parameter space location are presented for the six GUT models. In Section 
4, we determine the prospects for measuring the branching fractions for various 
Higgs boson decays, including those for specific supersymmetric (SUSY) sparticle 
pairs. The ability to discriminate between different GUT models and to determine 
the parameter space location within the correct GUT model on the basis of Higgs 
decays is delineated. Section 5 summarizes our results and conclusions. 

2 The GUT Models, Masses and Higgs Decays 

In the simplest GUT treatments of the MSSM, soft supersymmetry breaking 
at the GUT scale is specified by three universal parameters: 

• mo: the universal soft scalar mass; 

• mi/2: the universal soft gaugino mass; 

• Aq-. the universal soft Yukawa coefficient. 

The absolute value of /i (the Higgs mixing parameter) is determined by requiring 
that radiative EWSB gives the exact value of for the experimentally measured 
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value of mt] however, the sign of /i remains undetermined. Thus, the remaining 
parameters required to completely fix the model are 

• tan/3: the vacuum expectation value ratio; and 

• sign(/i). 

We remind the reader that a universal gaugino mass at the GUT scale implies that 
M3 : M2 : Ml ~ 3 : 1 : 1/2 at scale ~ rriz- For models of this class one also finds 
that \fi\ ^ Mi^2- These two facts imply that the Xi is mainly bino, while X2 
xt are mainly wino, with heavier gauginos being mainly higgsino 0. The running 
gluino mass m~{m~) is roughly three times as large as m~o ~ m~+ which in turn 
is of order twice as large as rn-o^. (The pole gluino mass is generally substantially 
larger than m~{m~) when squark masses are large.) 

We will consider three representative GUT scenarios characterized by increas- 
ingly large values of mo relative to mi/2 (which translates into increasingly large 
slepton masses as compared to m-o, m~o, and ^^+): 

• "No-Scale" (NS) @: = mo = 0; 

• "Dilaton" (D) §: mi/2 = -Aq = VSttiq; 

• "Heavy-Scalar" (HS): mo = mi/2, Aq = 0. 

Within any one of these three scenarios, the model is completely specified by values 
for mi/2, tan/5 and sign(/i). We will present results in the (mi/2,tan/3) parameter 
space for a given sign(/i) and a given choice of scenario. Our notation will be NS~ 
for the No-Scale scenario with sign(/i) < 0, and so forth. 

In Figures |I|, H and|^ we display the allowed (mi/2,tan/5) parameter space for 
the NS, D and HS scenarios, respectively. The boundaries of the allowed parameter 
space are fixed by experimental and theoretical constraints as follows: 

• The left-hand boundary at low mi/2 derives from requiring that Z ^SUSY 
decays not violate LEPl limits. 

• The low-tan/? boundary is obtained by requiring that the t-quark Yukawa 
coupling remain perturbative in evolving from scale m^ to the GUT scale. 

• In the NS scenario, the allowed parameter space is finite by virtue of two 
competing requirements. First, there is an upper bound on tan as a function 
of mi/2 obtained by requiring that the LSP (always the x? iri the allowed 
region) not be charged {i.e. we require m~ > ''^jo)-[| Second, there is the 
lower bound on tan /3 required by t-quark Yukawa perturbativity. One finds 
that for large enough mi/2 the upper bound drops below the lower bound. 

^This bound is especially strong in the NS scenario due to the fact that mo — implies very 
modest masses for the sleptons, in particular the ri, at scale mz- 
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The upper bound on tan as a function of mi/2 in the D scenario comes from 
demanding that the LSP not be charged. 

In the HS scenario, the upper bound on tan/3 arises by requiring that the 
SM-hke hght Higgs mass he above the current hmit of mi^a > 63 GeV. (In 
the HS scenario, for fixed '^a" becomes smaller and smaller as tan/3 

increases until eventually it approaches zero forcing m/jO to decline rapidly. 
In other scenarios, with lighter scalar masses and hence sleptons, the LSP 
becomes charged before tan/3 becomes so large that m^o starts declining 
rapidly.) 

In the D and HS scenarios, there is no upper bound on mi/2 unless cosmo- 
logical constraints are imposed. High mi/2 values (roughly, mi/2 ^ 500 GeV 
are, however, disfavored by naturalness considerations. 



Before proceeding, we provide a few technical notes. First, we note that the 
evolution equations must be implemented very carefully when considering very 
large m^o values. In order to avoid instabilities^ deriving from unnaturally large 
(and hence unreliable) one-loop corrections (for going from running masses to pole 
masses), we found it necessary to terminate evolution for soft masses at scales of 
order the associated final physical squark, slepton and heavy Higgs masses. In this 
way, the one-loop corrections are kept small and the physical masses obtained are 
reliable. The evolution program we employed is based on one developed by C. H. 
Chen |T^. Results at low mass scales were checked against results obtained using 



the programs developed for the work of Refs. and Once the appropriate 
low-energy parameters were determined from the evolution, we then employed 



ISASUSY [|rT| to obtain the branching ratios for the Higgs boson and subsequent 
chain decays. The ISASUSY results were cross-checked with our own programs. 
The decay results were then combined with Higgs boson pair production rates to 
determine rates for specific classes of final states. 



2.1 Sparticle and Higgs Masses 

Also displayed in Figs. |l], ||, and |^ are contours of constant rn-o, rn~+, mj^, 
and m^o. These reveal the importance of detecting the heavy Higgs bosons and 
measuring their masses accurately. The masses of the inos and the sleptons will 
presumably be measured quite accurately, and the figures show that they will 
determine in large measure the values of mi/2 and mg. But the rather vertical 
nature of the m-o, '"^^±5 and mj^ contours implies that tan/3 is likely to be poorly 
determined from these masses alone. Fortunately, the m^o contours are not nearly 
so vertical, implying that a measurement of m^o can be combined with the mi/2 
determination from the ino masses to fix a value of tan/?. The accuracy of this 

^Such instabilities are found, for example, in ISASUGRA 
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No Scale Scenario Mass Contours (GeV) 




Figure 1: We show the {mi/2,'tsin f3) parameter space regions (bold outer 
perimeter) within which we find a consistent EWSB solution for the No- 
Scale model. Contours of constant mass are shown within the allowed 
region for the Xi, xf^ ^^id Ir. Results for both signs of n are shown. 
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Figure 2: We show the (mi/2,tan P) parameter space regions (bold outer 
perimeter) within which we find a consistent EWSB solution for the Dilaton 
model. Contours of constant mass are shown within the allowed region for 
the Xi, xfy and Ir. Results for both signs of fj, are shown. 
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Heavy Scalar Scenario Mass Contours (GeV) 
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Figure 3: We show the {mii2-,tan (3) parameter space regions (bold outer 
perimeter) within which we find a consistent EWSB solution for the Heavy- 
Scalar model. Contours of constant mass are shown within the allowed 
region for the xf , and Ir. Results for both signs of ji are shown. 
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determination depends upon the accuracy with which m^o (and mno, TnH±) can 
be measured. For discovery in the — >■ bb decay mode (as possible for almost 
all model parameter choices at full luminosity, see later), this accuracy is fixed by 
the bb mass resolution. At an e~^e~ collider, a resolution of itAMbf, ~ ±10 GeV 
is probably attainable. For a large number, N, of events, m^o can be fixed to 
a value of order AMi,b / VN, which for = 20 (our minimal discovery criterion) 
would imply Am^o ~ 2 — 3 GeV. Examination of the figures shows that such mass 
uncertainty will lead to a rather precise tan (3 determination within a given GUT 
model, except at low ttiao and high tan /? in the NS case. 

2.2 Higgs Decays 

Let us now turn to the decays of the heavy Higgs bosons of the MSSM. As 
already noted, our ultimate goal is to use these to confirm/re-enforce the correct- 
ness of both the model and the parameter choices within the model that has been 
singled out by the mass measurements. The most important common feature of 
the GUT models we consider is that squarks are always sufficiently heavy that 
decays of Higgs bosons to squark pairs are not kinematically allowed. This is true 
even for the NS boundary conditions with rriQ = 0, in which the large squark 
masses derive from the substantial evolution of the colored soft-scalar masses to 
positive values as the scale decreases from Mjj towards mz- In order that the 
squarks be light enough for squark pairs to appear in Higgs decays, substantial 
breaking of the universality of soft-SUSY-breaking scalar masses at the GUT scale 
is required. For example, light sbottom and stop squarks can be consistent with 
radiative EWSB via evolution if the Higgs soft scalar masses are much larger than 
the squark (in particular, stop and sbottom) soft scalar masses at Mu. In this 
case, if°,y4° — > titi,bibi and —>■ tibi pair channels would dilute the SM decay 
modes of the Higgs to a much greater extent than do the ino and slepton decays in 
the models discussed here. Strategies for detecting and studying H^A^ and H^H^ 
pair production would have to be reconsidered. In any case, there would be no 
difficulty in distinguishing models with light stop and/or sbottom squarks from 
the NS, D and HS models considered here. 

In fact, the three models we consider are rather similar to one another in most 
respects. Thus, they provide a good testing ground for assessing the extent to 
which we can distinguish between models by using experimental information from 
the Higgs sector. We shall see that Higgs branching ratios depend substantially on 
the particular model choice and on the precise location in parameter space within 
a given model. Figures and ^ illustrate the dependence of Higgs branching 

fractions upon parameter space location for the yU < Dilaton (D~) scenario. In 
these figures, we give contours of constant branching fractions for A^ and 
decays. The decay channels 66, tt, XiXi) X1X2; the sum over all SUSY 
decay channels, are considered for both the and A^. In addition, we show the 
h^hP and VV (summed over all V types) branching fractions for the H^. (The VV 
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Dllaton (/2= 



— ) Branching Ratio Contours 




Figure 4: a) We show contours within the {1711/2, tan (3) parameter space of 
constant branching fraction for the bb, ti, xtxi ^ hPhP, and X1X2 

decay channels. Results are for the D~ scenario. 
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Figure 4: b) We show contours within the (mi/2,tan/3) parameter space 
of constant branching fraction for the — > bb, tt, xtxij X1X2 decay 
channels, as well as for — SUSY and — SUSY, summed over all 
SUSY channels. Results are for the D~ scenario. 



12 



Dllaton (/2= 



— ) Branching Ratio Contours 




Figure 4: c) We show contours within the (mx/2)tan/3) parameter space 
of constant branching fraction for the tb, W^h^, i^u, t'^u, cs, and 

SUSY, summed over ah SUSY channels. Results are for the D~ scenario. 
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branching fraction for the is very tiny.) The A° Zh^ branching fraction is 
small, but, as we shall see, measurable in some regions of parameter space. For 
the H'^ we display branching fraction contours for tb, W~^h^, iiu, t'^u, cs, and the 
sum over all SUSY decays. [B{A° Zh^) is similar to B{H+ W+hP)] Several 
important features of these plots deserve emphasis. 

• For the and A^, the net branching fraction for SUSY decays dechnes 
rapidly with increasing tan (5 due to the enhancement of the bb coupling and, 
hence, increasing relative importance of bb decays. 

• SUSY decays of the and A^ are also small when mHO,mAO > 2mt, with 
the relative branching fraction 5(SUSY) / B{tt) saturating to a constant value 
below 0.1 for large m^o (equivalently large mi/2) at fixed tan/3. 

• For mHOjmAO > 2mt, the ratio of bb to tt branching fractions rises very 
rapidly as tan /3 increases. 

• The SUSY decay branching fraction of the is relatively independent of 
tan /3 for lower mi/2 values. 

• B{H^ —>■ W^h^) [as well as B{A^ — >• Zh^)] is only signficant when tan/3 and 
mi/2 are both small. 

• B{H^ h^h^) is significant for a larger range of modest tan/3 and mi/2 
values than the former two branching fractions. 

• B{H^ T^v) remains significant (> 0.1) for a range of tan/3 values that 
becomes increasingly large as mi/2 increases. 

These figures show that a measurement of several ratios of branching fractions {e.g. 
SUSY/66 for the H^.A^ and SUSY/t6 for the H^) would determine the values of 
tan/3 and mi/2- Branching ratios in the other five scenarios display a more or less 
similar pattern to that found in the D~ case, although the numerical values at any 
given (mi/2,tan/3) location can differ substantially. For any given GUT scenario, 
definite predictions for all other experimental observables are then possible and 
could be checked for consistency with observations. In particular, the predicted 
Higgs, ncutralino, and chargino masses should agree with the measured values if 
the GUT scenario is the correct one. 

3 Discovering the i/^, and 

In this section, we determine the luminosity required in order that discovery 
of H^A^ and H^H~ be guaranteed over essentially all of the allowed parameter 
space of the three scenarios. For the models considered in this paper, we find that 
discovery is always easiest by employing final states in which neither of the Higgs 
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bosons of the pair decays to a final state containing SUSY particles. The final state 
configurations we employ for discovery are listed below, along with techniques for 
isolating them from backgrounds. 

• I) H'^A^ — > 46: We demand observation of four jets which separate into two 
nearly equal mass two-jet pairs. Event rates for this mode [labelled by N{4b)] 
include a factor of B{H^ bb)B{A^ bb). 

• II) H°A° with i/o ^ h^h° 46 and A^ X: it would be sufficient to 
observe the two /i°'s by demanding two jet pairs that reconstruct to the 
known m/jo recoiling against a reconstructed (from incoming energy and net 
h'^h'^ pair four momentum) 'missing' mass that is the same as the h^h^ pair 
mass. Event rates for this mode [labelled by N[hh)] include a factor of 
5(1^0 ^/i0/i0)[5(/i0^66)]2. 

• III) H'^A'^ 4t: We can simply demand > 10 visible (and moderately ener- 
getic/separated) leptons/jets. The predicted rate for such states on the basis 
of QCD (including it production) is quite small. Because of inefficiencies 
associated with combinatorics, we would not require direct reconstruction of 
the 1^'s or f s (implying that we would also not be able to require roughly 
equal Higgs boson masses). Event rates for this mode [labelled by N{4t)] 
include the effective branching ratio for — ^> 10 visible leptons/jets, 
given by B(H° tt)B{A° tt)B{tttt ^> 10 visible). 

• IV) H^H^ 2t2b: We insist on 8 jets or 1 lepton plus 6 jets (in particular, 
fewer than 10 visible leptons/jets so as to discriminate from the above 4t final 
states) and possibly require that one W and the associated t be reconstructed. 
Event rates for this mode [labelled by N{tb)] include a factor of [B{H^ 
tb)]'^ {2B{t 2jb)B{t e+ub) + [B{t 2j6)]2}. 

There will also be an overall efficiency factor for detector coverage and for ex- 
perimentally isolating and detecting these modes. This will be incorporated in 
our yearly event rate estimates by reducing the total luminosity available (pre- 
sumed to be L = 200 fb^^ per year at ^/s = 1 TeV and L = 1000 fb^^ per year 
at ^/s = 4 TeV) by an overall efficiency factor of 40% (to L^s = 80 fb^^ and 
Leff = 400 fb~^, respectively). We have not performed a detailed simulation, but 
beheve that that such an efficiency is not unreasonable given the fact that back- 
grounds are relatively small for the above outlined signatures. In particular, since 
all the final states contain at least four b jets, we can require one or two 6-tags 
(in order to eliminate any residual QCD background) without incurring significant 
penalty, given that the vertex-tagger should have efficiency of 60% or better for 
any single 6-jet within its acceptance. (Tagging of 6-jets would be desirable for 
cleanly separating H^A^ from H^H~ final states. In the absence of any 6-tagging 
there would be a small number of H^H~ — > 4j (with j = c, s) events that would 
combine with the H^A^ — > 46 final states to the extent that mH+ ~ tuho ~ m^o.) 
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After including branching fractions and the 40% efficiency, something hke 20 events 
should be adequate for detection. In our graphs, we will display 20, 50 and 200 
event contours. 

If the H^A^ —>■ 4t mode is dominant, we will wish to reconstruct the mass of 
either the or the from the 4j26 decay mode of one of the tt pairs. This 
will be important both as a means for measuring the mass and also as a means 
for triggering on H^A^ pair production using just one of the two members of the 
pair (see Section 4). There will be a further efficiency factor (on top of the above 
overall efficiency factor) for isolating the relevant events and then reconstructing 
the mass of the A° or H^. We estimate this additional efficiency factor be of order 
25% each for the A^ and H^. This is the result that would be obtained from 
[B{t 2j6)]^e, with e = 0.55 for combinatoric and other problems. The low net 
efficiency (~ 0.2 = 2x ~ 0.25 x 0.4) for events in which either mj^o or m^o could 
be fully reconstructed imphes that an accurate determination of m//o ~ m^o would 
require several years of running if H^A^ At is the dominant final state. 

There are several reasons why non-SUSY final states are best for discovery: 

• As illustrated in Figs. |a-c, branching fractions for SM decays, e.g. A^, —>■ 
bb or tt and tb, do not fall much below 0.1; 

• Unlike the bb channel, mass reconstruction in SUSY modes is not possible 
(due to missing energy). 

• Particle multiplicities in the At and 2t2b final states are sufficiently large to 
be very distinctive and free of background, unlike many of the final states 
associated with SUSY decays. 

In Figures ^, H, and ^ (for the NS, D and HS scenarios, respectively) we give 
the 20, 50 and 200 event contours in the (mi/2,tan/?) parameter plane for H'^A^ 
discovery modes I, II and III and H^H^ discovery mode IV at ^/s = 1 TeV. We 
assume L = 200 fb~^ and e = 40% efficiency, i.e. Lcs = Le = 80 fb~^. Results 
are displayed for both signs of fi. Also shown are the boundaries defined by the 
kinematically accessible rriHO + rriAo < y/s or 2mj/± < y/s portion of the allowed 
parameter space (bold solid lines). In comparing scenarios, it will be important to 
note that the NS scenario plots have greatly expanded axis scales relative to plots 
for the D and HS scenarios. 

As noted earlier, 20 events is likely to be adequate for discovery; the 50 event 
contour at Leg = 80 fb~^ would probably allow discovery at Leg = 32 fb~^ and the 
200 event contour would allow discovery at Leg = 8 fb~^. These figures show that 
for all three GUT scenarios at least 20 H'^A^ events are present in one or more of the 
modes I-III throughout almost the entire kinematically accessible portion of the 
allowed (mi/2,tan/5) parameter space. If the 50 event contours are appropriate 
(because Lgfr is a factor of 2.5 smaller) then one begins to see some, but not 
enormous, sections of parameter space such that H^A^ detection would not be 
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possible. If efficiencies and integrated luminosity were in combination a factor 
of 10 worse than anticipated, the 200 event contours might apply; they indicate 
that H^A^ detection would then be possible only in the part parameter space 
characterized by small values of mi/2 and large values of tan/5. 

In these same figures, the 20, 50 and 200 event contours for the H^H~ discovery 
mode IV show that 20 events are found for all of the constraint and kinematically 
allowed parameter space except a small wedge at small mi/2 values. The 200 event 
contours (equivalent to 20 events at L^s = 5 fb~^) cover nearly as a large section of 
parameter space. Thus, even if efficiency and luminosity are in combination a factor 
of 10 worse than anticipated, H^H~ discovery after just one year of running would 
be possible over the bulk of parameter space. The somewhat better guarantees for 
the H^H~ mode as compared to the H^A^ mode derive simply from the larger 
H^H~ cross section which is roughly a factor of 3 larger than that for H^A^. 

This same analysis can be repeated for 4TeV and L = 1000 fb~ (imply- 
ing L = 400 fb~^ for e = 0.4 efficiency) with very similar results. The kine- 
matic range is much greater, allowing H^A^ and H'^H~ production out to masses 
m^o ~ rriHO ~ rnH± < 2 TeV. (The limited NS parameter space implies that such 
energies are not needed were this the correct GUT scenario). If 20 events are ade- 
quate (and they would certainly be rather spectacular events at high Higgs boson 
mass) then both H^A^ and H^H^ detection would be possible for nearly all of 
the constraint/kinematically allowed parameter space for all three GUT scenarios. 
Dimunition in coverage due to poorer efficiency or lower luminosity follows much 
the same pattern as described for the a/s = 1 TeV, L = 200 fb~^ case. To illus- 
trate, we present the 20, 50 and 200 event contours for H^A^ (final states I-III) 
and H^H~ (final state IV) in the D scenario. Fig. ^. 

4 Measuring Ratios of Branching Fractions and 
Discriminating Between Models 

In this section we discuss the prospects for measuring the relative size of the 
various branching fractions for different decay modes of a given Higgs boson and 
for using such measurements to pin down the GUT model and parameter choices 
within a given GUT model. Additional information is contained in the absolute 
rates for different types of final states. However, it is likely that greater uncertainty 
will be associated with absolute rates than with ratios of rates, since some types 
of efficiencies will cancel out of the ratios. 

The key to determining the relative magnitude of the branching fractions for 
different final state decays is to first identify and mass-reconstruct ('tag') one of 
the Higgs bosons in the H^A^ or H^H^ pair final state, and then compare the 
relative rates for different types of decays of the second Higgs boson. Identification 
and mass-reconstruction of the first Higgs boson requires using one of its fully 
reconstructable final states. As additional verification that the event corresponds 
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Figure 5: We show NS model 20, 50 and 200 event contours within the kine- 
matically accessible portion of the allowed (mi/2, tan /?) parameter space for 
H^A^ discovery modes I, II, III and H^H^ discovery mode IV, assuming 
Leff = 80 fb-\ 
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Figure 6: We show D model 20, 50 and 200 event contours within the kine- 
matically accessible portion of the allowed (mi/2, tan /?) parameter space for 
H^A^ discovery modes I, II, III and H^H^ discovery mode IV, assuming 
Leff = 80 fb-\ 
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Figure 7: We show HS model 20, 50 and 200 event contours within the kine- 
matically accessible portion of the allowed (mi/2, tan /?) parameter space for 
H^A^ discovery modes I, II, III and H^H^ discovery mode IV, assuming 
Leff = 80 fb-\ 
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Dlllaton Scenario Discovery Contours 
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Figure 8: We show D model 20, 50 and 200 event contours within the kine- 
matically accessible portion of the allowed (mi/2, tan /?) parameter space for 
discovery modes I, II, III and H^H^ discovery mode IV, assuming 
Leff = 400 fb"^ at ^/s = 4TeV. 
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to Higgs pair production, we would require that the missing mass (as computed 
using the incoming center-of-mass four-momentum and the four-momentum of the 
reconstructed Higgs) be roughly equal to the mass of the identified Higgs. For 
identification and mass-reconstruction of the first Higgs boson, we employ: 

• H^A^ with 2b or A° ^ 26; 

• H^A^ with H'^ — » 2t or A° — > 2t — note that, unlike the 4t discovery channel, 
reconstruction of the 2t mass will be necessary, and will be accompanied 
by an extra efficiency penalty relative to 2b or A^ 2b tagging of 
ett/bb = [B(t 2j6)]^e ~ 0.25 (for e = 0.55), as discussed earlier; ^ 

• H^A^ with h^h^ 46; 

• H^H^ with ^ tb W2b 2j2b, or the reverse — tb mass reconstruc- 
tion will be necessary. 

In the case of H^A^ pair production, in determining that the second (non-tagged) 
member of the pair decays to tt, we will again demand full ti reconstruction, and 
we will apply the extra eu/bt efficiency penalty relative to bb decay. This might 
be somewhat too conservative an approach, but does simplify our analysis since 
the event rates of interest involving tt + bb decays will then be proportional to the 
effective branching fractions 

5eff (if°, A'^bb + tt) = B{H\ AO ^ bb) + eu/bbB{H\ A° ^ tt) . (1) 

Because m^io ~ tuho over much of parameter space, we will presume that it 
is not possible to separate the A^ and from one another. We also stick to 
our simplifying assumption that the overall efficiency, e, associated with detector 
coverage, 6-tagging and so forth does not depend upon the final state, except 
that in the case of it decay we include an extra e in eu/bbi as discussed above 
and as incorporated through B^e defined in Eq. ([l|). With these assumptions, the 
following ratios of branching fractions can be extracted directly from experimental 
observations using the measured values of B{h^ bb) and B{t 2jb). 

B{H^ ^ SUSY)5cff(A0 ^bb + tt) + B{A^ SUSY) B,s{ H^ bb + tt) 



5eff(i/o ^bb + tt)B,s{A^ ^bb + tt) 
B{H^ -> tt)B{A° -> bb) +B{A° ^ tt)B{H^ bb) 
B{H^ bb)B{A^ bb) 
B{H^ h^h^)B{A^ bb) 



(3) 



(4) 



B{H^ bb)B{A^ bb) 

These details for the ti final state are only relevant for BeS defined in Eq. (j^) and the ratios 
of Eqs. (H) and (|^) below, and then only when B{H^ , ti) are relatively large. 
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(9) 



As a shorthand, we will employ the notations 

/2B{H'^,A^ ^ SUSY)\ /2B{H^,A^ ^tt)\ 
\ B,siH<^,A<^ ^bb,tt) I ' \ B{H^, AO bb) I ' 

for the ratios of Eqs. (|]) and @, respectively. The ratios of Eqs. (H)-® reduce to 

2B{H+ SUSY, r+z/, W+h^) 



B{H+ tb) 



(10) 



respectively. We retain both hb and tt final states in Eq. (^, using the combination 
defined in Sefr, in order that we may assess the importance of SUSY decays both in 
regions where hb decays of the if", A^ are dominant and in regions where ti decays 
are important. 

In estimating the accuracy with which these ratios can be measured experimen- 
tally, it is important to keep track of the actual final state in which the observation 
occurs and the effective efficiency for observing that final state. We make this 
explicit below. 

• The event rate for the numerator of Eq. (j^) is obtained by multiplying the 
rate for H^A^ pair production by e (the overall efficiency factor) times the 
indicated sum of branching ratio products: [B{H^ SUSY)5cfr(^'' — >■ 66 + 
tt) + B{A^ SUSY)5eff(if° ^bb + tt)]. 

• The numerator of Eq. (|) must be measured in the final state in which both 
/i°'s decay to hb. Thus, the event rate associated with determining the numer- 
ator is obtained by multiplying the H^A^ pair production rate by a factor of 
[B{h^ bb)f times the overall efficiency e times B{H^ h^h^)B{A^ bb). 

• The event rate associated with measuring the numerator of Eq. (^ is obtained 
using a factor of eeu/bb = ^e[B{t 2jb)]^ times B{H^ tt)B{A° bb) + 
B{A° tt)B{H° bb). 

• The event rate for the numerator of Eq. (^) is computed using the factor 
eB{h^ bb)B{A° Zh^)B{H° bb). This implicitly assumes that we 
can sum over all Z decays, as would be possible since the Z mass can be 
reconstructed from the cm. ^/s value and the momenta of the four 6's. 
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• The event rate for the numerator of Eq. (|^) is obtained by multiplying the 

H+H- event rate by the factor eB{t 2jb)[B{H+ SUSY)B{H- 
bt) + B{H- ^ SUSY) B{H+ ^tb)]. 

• The event rate for the numerator of Eq. (|^ is computed by multiplying 
the pair rate by eB{t 2jb)[B{H+ t+u)B{H- bt) + B{H- 
T-v)B{H+ tb)]. 

• The rate for the numerator of Eq. (|^) is computed using the factor eB{hP — »• 
bb)B{t 2jb)[B{H+ h^W+)B{H' hi) + B{H- h^W-)B{H+ 
tb)]. 

• The factors for the denominators are obtained by multiplying the indicated 
branching ratio product by e in the case of the neutral Higgs ratios, and by 
e[B{t 2jb)Y in the case of charged Higgs ratios. 

• When dividing the SUSY collection of final states (as simply identified by 
missing energy) into subcategories of a certain number of leptons and/or 
jets, the full set of appropriate branching ratios are included in all the chain 
decays leading to the specified final state. 

As noted, the overall factor of e common to all rates is incorporated by reduc- 
ing the full luminosity to the effective luminosity Leg. Rates for the standard 
Leff = 80 fb~^ are thus obtained by computing the pair production cross section, 
multiplying by L^s and then including all the above factors after removing the 
overall multiplicative e contained in each. The bottom line is that even though we 
plot the ratios listed, the statistical errors we shall discuss will be based on the 
actual number of events as obtained according to the above-outlined procedures. 

The utility of the above ratios derives from the following general features. The 
1st ratio is primarily a function of tan/3. The 2nd provides an almost direct 
determination of tan/5 since is roughly proportional to cot^/? in the MSSM. 
The ratios of Eqs. (^) and (0) both exhibit substantial and rather orthogonal 
variation as a function of tan/5 and mi/2. The ratio of Eq. is proportional 
to the relative strength of the —>■ h^h^ trilinear coupling as compared to the 
bb coupling. This could be the first direct probe of Higgs trilinear couplings. 
The ratios of Eqs. (||) and (|]) would probe the very interesting Higgs-Higgs-vector- 
boson couplings. These features will be illustrated shortly. 

4.1 Resolving Ambiguities in Identifying Different Final 
States 

Since all SUSY final states will contain substantial missing energy, the ambi- 
guities in separating SUSY decays from others are limited. We discuss below the 
procedures for removing the only ambiguities that appear to be of importance. 
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(A) A potential ambiguity arises in + t final states of the to which the 
SUSY and x'^x^ decay modes and the SM — >■ r"*"!/ — >■ £~^3i' decay modes 
all contribute. The —>■ t^Pt —>■ i'^Su decay can be identified using kinematic 
constraints. Consider the cm. system of the decaying if+ (as determined using 
incoming beam information and the tagged H~ four-momentum). To the extent 
that rrir can be neglected and, therefore, the r decays coUinearly to f"'"2z/, all of 
which move opposite the primary i/^-, one must have E = \^t\, where E is the 
energy of the observed i. SUSY events of any type will normally violate this 
constraint. In what follows, the i'^V and x^X° decays are both included in the 
overall SUSY decay rate of the 

(B) In or decay, t'^t~ decays contribute to the same £'^£~ + final 
states to which the SUSY modes and x~^x~ contribute. The procedure for 
eliminating the r+r" decay is analogous to that discussed in (A) for removing 

— > t'^i't decays. We again note that, for most events, the r mass can be 
neglected relative to its momentum. In the (known) rest frame of the Higgs, the 
coUinear approximation implies that the and £~ and their associated neutrinos 
travel in essentially the same directions as the parent r"*" and r~, respectively. As 
a result, such events must have \E^ — E_\ = \I^t\, where E± are the observed 
energies of the in the Higgs rest frame. The very non-coUinear SUSY modes 
would generally be far from approximately satisfying this constraint. Kinematic 
constraints do not allow an event-by-event separation of the two SUSY modes, 
£~^£~ and x'^x~^ in the + final state. These are lumped together as part 
of the overall SUSY decay branching fraction. 

(C) Events in which the unreconstructed Higgs boson/decay is orA^ ^ ti ^ 
£v2]2b or ^ tb ^ £vbb can be eliminated by using the incoming beam en- 
ergy/momentum 4- vector, subtracting the momenta of all visible final state lep- 
tons and jets, and computing the invariant mass of the resulting difference 4- vector. 
This would belong to the v in the above cases. A cut requiring a substantial value 
would eliminate the above final states and be highly efficient in retaining true 
SUSY decays. For parameters such that the rates for single neutrino events (as 
defined by the above procedure and requiring a small value for the difference 4- 
vector mass) are significant, we shall find that the H^,A^ — > tt and — > ib 
branching fractions can be directly measured with reasonable accuracy (using all- 
jet modes). The predicted single neutrino rate could then be compared to that 
observed as a further check. Events where the unreconstructed Higgs/decay is 

or A^ tt 2£2v2b cannot be eliminated by the above technique. However, 
the branching fractions, B{H^,A^ tt), measured in all-jet final states can be 
employed to make an appropriate correction. The single neutrino rates, as defined 
above, may allow a double-check of the all-jet final state determinations of the it 
branching fractions. 

(D) Other ambiguities include events in which the second Higgs boson/decay is: 
A^ — >• ZhP — >• Zr^r~ , where the Z decays invisibly or to r+r^; — >• W^hP — >• 
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W^T'^T~ where the decays leptonically; and — > h^h^ — > t~^t~t^t^ . The 
common characteristic of all these is the presence of missing energy from r, W 
and/or Z decays that is due to more than a single neutrino and that makes it 
impossible to either directly or indirectly reconstruct the mass of the h^, W and/or 
Z. However, the event rates for these processes are so low that they can be included 
in SUSY decays without any visible alteration of the effective SUSY branching 
fraction. Further, whenever the — > h^h^, —>■ Zh^ or — »• W^h^ decays 
are significant, we shall see that at least a rough measurement of the corresponding 
branching ratio will be possible in all-jet modes. Given the known Z vV and 
/i" T^T^ branching fractions, a correction could then be made using a Monte 
Carlo simulation. 

4.2 Ratio Contours, Error Estimates and Model Discrimi- 
nation 

In order to determine how well we can measure the ratios of Eqs. (H), (j^), @, 
(§)> (@)> (0)' ®' have proceeded as follows. For each of the six scenarios 
(D^, D+, . . . ) and for a given (mi/2,tan/3) choice within the allowed parameter 
space of a given scenario, we first compute the expected number of events available 
for determining the numerator or denominator of each ratio. The ingredients (such 
as branching ratios and efficiencies) in the event number computations for each 
channel were given earlier. The expected number of events in the numerator or 
denominator is taken as the mean value in determining a Poisson distribution 
for that event number; if the mean number of events is > 30, then we use a 
Gaussian approximation to the distribution. From the event number distributions 
we compute the probability for the numerator and denominator of each ratio to take 
on given values. (We fluctuate the event numbers and then correct for branching 
fractions and efficiencies.) The probability of the resulting value for the ratio 
is then simply the product of these probabilities. The probabilities for different 
combinations that yield the same value for the ratio are summed. In this way, we 
obtain a probability for every possible value of the ratio. These probabilities are 
re-ordered so as to form a distribution. The lower (upper) limit for the ratio at 
this (mi/2,tan/9) value is then found by adding up the probabilities, starting from 
zero, until the sum of is 15.9% (84.1%). In other words, the confidence level that 
the true value of the ratio is higher (lower) than the lower (upper) limit is 84.1%. 
These would be the ila upper/lower limits for the ratio in the limit where the 
distribution of the ratio is normal. 

In computing the number of events available for determining the numerator 
or denominator (or one of the independent contributions thereto) we include only 
fully reconstructable final states for the tagged Higgs boson. The branching ratios 
and efficiency factors were detailed below Eq. ([I0|). We presume Lgfj = 80 fb~^ at 
1 TeV (Leff = 400 fb^^ at 4TeV). The efficiency factor, e, included in L^^ should 



26 



reflect efficiencies associated with identifying a particular type of event in such a 
way as to eliminate backgrounds, e.g. via 6-tagging, cuts on I^t, and so forth; 
the e appropriate to the current situation where one of the Higgs must be clearly 
'tagged' (as defined earlier) will probably be smaller than that appropriate to sim- 
ply discovering a signal, given the need to clearly separate different types of final 
states from one another. Thus, the above Lgg = 80 fb~^ (400 fb~^) values probably 
would only be achieved after several years of running. We re-emphasize that an 
implicit approximation to our approach is that L^s is the same for all the obser- 
vationally/statistically independent final states that appear in the numerator and 
denominator of a given ratio. [Aside from our special e = 0.55 correction for tt re- 
construction, the only explicitly channel-dependent factors that have been included 
are the relevant branching fractions, as detailed below Eq. (|10]).] Presumably, this 
will not be true in practice, but it is at least a reasonable ffist approximation. Full 
detector specification and simulation would be necessary to do better. 



In Figs, pi ro, O, O, 13, 14, and O we plot contours of constant values for 



the ratios of Eqs. (0), (§), (|), (§), (§), (0), and (§ within the ^/s = 1 TeV 
constraint/kinematically allowed (mi/2,tan/5) parameter space. Associated with 
each such contour, we give two additional contours showing how much the tan/3 
value at a given (known) value of mi/2 would have to change in order to reproduce 
the values obtained for deviations in the ratio at the ±lcr statistical level. [As 
previously explained, ±lcr is our short hand phrase for deviations such that the 
ratio has 84.1% probability of being lower (higher) than the upper (lower) limit.] 
We do not consider errors when there are fewer than 4 events that can be used to 
determine the numerator for one of these ratios. The 4-event contours are indicated 
on the figures. 

Consider ffist the relative SUSY branching ratio contours of Eqs. (0) and (^) 
displayed in Figs. ^ and |13|, respectively. For most points in parameter space, 
a simultaneous measurement of the two ratios will determine a fairly small and 
unique region in the parameter space of a given model that is simultaneously 
consistent with both measurements at the la level. 

If tan/3 is not large, then measuring B{H~^ t^v) / B{H^ tb) via the 
ratio of Eq. (|^ can provide a second determination of tan p. The dependence of 
this ratio on tan/? for a selection of mH+ values is illustrated in Fig. |1^. There, 
the ratio is computed at tree level. We see that the ratio depends sensitively on 
tan/3 at fixed mH+ for tan/3 < 6. For such tan/3 values, measurement of the ratio 
provides an excellent tan /? determination. However, when tan (3 is large the r^v /tb 
ratio becomes independent of tan/? and sensitivity is lost. Note also that the ratio 
becomes independent of mH+ when mH+ is large. Thus, when both mH+ and tan j3 
are large, this ratio will provide little information regarding location in parameter 
space. 

Contours of constant B{H^ — > t^v) / B{H^ tb) in (mi/2,tan/3) parameter 
space are displayed in Fig. |l^. It is also useful to plot these same contours in 
(m//+,tan/3) parameter space, as done in Fig. O. In both figures, one observes 
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2<B(H°, A°^SUSY)/B,„(H°, A°^bb or tt)> 



e^e'^H^A", E,^=1 TeV, U=80 fb" 
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Figure 9: We plot contours, along which the ratio of Eq. (|2|) has a given 
constant value, within the constraint/kinematically allowed (mi/2) tan /?) 
parameter space of the D~, D"*", NS~, NS"*", HS~, and HS^ models. Results 
are shown for the same three central values for all models. For each central 
value, three lines are drawn. The central line is for the central value. The 
other two lines are contours for which the ratio deviates by itlcr statistical 
error from the central value. Bold lines indicate the boundary beyond 
which fewer than 4 events are found in the final states used to measure the 
numerator of the ratio. 
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Figure 12: As in Fig. |9|, but for the ratio of Eq. (^). 
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2*B(H^^SUSY)/B(H*^tb) 




Figure 13: As in Fig. |9|, but for the ratio of Eq. (^). 
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2*B(H^ ^ tV)/B(H^ ^ tb) 




2*B(H*^h°W^)/B(H^^tb) 
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Figure 16: The ratio B{H^ — t^v)/B{H'^ — tb) computed at tree level 
for mt = 175 GeV and ruf, = 4 GeV as a function of tan/? for mfj+ = 200, 
300, 400, 600, and 1000 GeV. 
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2*B(H* ^ tV)/B(H" -> tb) 
eV^H^H", E,^=1 TeV. U=80 fb 
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Figure 17: Contours of the ratio B{H+ t+v)/B{H+ tb) and its 
associated ±la contours are plotted as a function of tan/? and for 
Leff = 80 fb"^ at ITeV. 
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2<B(H°, A°^tt)/B(H°, A°^bb)> 

H°A^ E,„=1 TeV, U„=80 fb=^ 
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Figure 18: Contour of the ratio 2{B(R^,J^ tt)/B{H^,A^ bb)) its 
associated ±la contours are plotted as a function of tan/3 and m^o for 
Leff = 80 fb"^ at 1 TeV. 
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a change from horizontal to vertical contours as one moves from low tan/5 and 
large mi/2 (equivalent to large mH+) to high tan/3 and small mi/2 (implying small 
mH+). The horizontal nature of the contours at large mi/2,'rnH+ and small tan/3 



can be understood from Fig. |I6[ As already briefly noted, this figure shows that 
when tan (3 is small, small changes in tan (3 yield large changes in the ratio, whereas 
there is little sensitivity to changes in mH+ at fixed tan/3 when mH+ is large. In 
contrast, for small mH+ Fig. shows that small changes in mH+ produce large 
changes in the ratio, whereas there is almost no sensitivity to tan/5 when tan/5 
is large. As a result the contours in Figs. |H| and are vertical at small mH+ 
when tan/3 is large. The wide separation between the central and ±la contours 
when mi/2 and tan/3 are both large is a reflection of the constancy of this ratio 
(as displayed in Fig. |T6|) when both tan /? and mH+ are large. Outside the region 
where tan/? and 771^+ are both large, the T^v/tb contours are roughly 'orthogonal' 
to those for the two SUSY ratios discussed earlier. 

In general, it is apparent that the contours for the ratios of Eqs. (^, (|^) and (0) 
in the (mi/2,tan/5) plane are all oriented rather differently. This means that, in 
combination, these three relative Higgs branching fractions provide a fairly power- 
ful check of the consistency of a given model, as well as a very definite determination 
of the value of tan/? that is required for a particular value of mi/2 in the model. 
We have already noted that mi/2 will be accurately determined in a given model 
by the neutralino and chargino masses, and that the measured m^o will generally 
provide a tan (3 determination. This determination of tan j3 from the masses and 
the value for tan (3 required for consistency with the above three ratios of branching 
fractions are usually not consistent with one another for an incorrect model choice. 

Additional discrimination power between the correct and an incorrect model 
choice is possible if we resolve the SUSY rates in Eqs. (0) and (H) into final states 
with a fixed number of leptons plus any number of jets (including 0) plus missing 
energy. Thus, instead of the single ratio of Eq. (^, where SUSY was defined to be 
the sum over all supersymmetric decay channels, it will prove useful to consider 
the three ratios obtained by dividing SUSY into the (i) [0^][> Oj], (ii) [U][> Oj] 
and (iii) [2^][> Oj] channels, where the [> Oj] notation indates that states with 
any number of jets (including 0) are summed over. Rates with [> 3£][> Oj] are 
negligible. Similarly, instead of the ratio of Eq. (|^) we will consider the two ratios 
obtained by separating SUSY into the channels (i) or (ii) defined above. Rates with 
[> 2£][> Oj] are negligible. All SUSY final states will have large missing energy. 
The five observable SUSY ratios so obtained are not very closely correlated, and 
thus are unlikely to be consistent with one another and with the t^v /tb ratio for 
any but the correct model choice. 

Still more discrimination power can be achieved via the other branching fraction 
ratios defined in Eqs. (H), (||), (|), and (§). For example, we see from Fig. |10 



that the tt/hh ratio is quite sensitive to tan/3. This is even clearer by displaying 
the contours in (m^o,tan/?) space. Fig. |18[ The h^h^/H^ bb, — > 



Zh^A^ bb and H+ W+h^/H+ tb ratios plotted in Figs. ITT], p] and 
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p!5| , respectively, are also sensitive to tan/?. However, even more interesting is 
their sensitivity to the sign of the /i parameter. All three ratios are much smaller 
for n > than for /i < [at a fixed (mi/2,tan/?) location]. These differences 
derive almost entirely from a large decrease in the h^h^, Zh^ and 

W~^h^ couplings, respectively, as the sign of fi is changed from + to — . (In 
the case of the —>■ h^h^ coupling, this decrease is largely due to the change of 
sign of a radiative correction to the vertex associated with top, bottom, stop and 
sbottom loops. In the —>■ Zh^ and — > W~^h^ cases, the large decrease is a 
tree- level effect.) Together, these three ratios will provide significant discrimination 
between scenarios with the opposite sign of /x. 



4.3 Quantitative Strategy for Estimating Model Discrimi- 
nation Power 

To determine the discrimination power achieved by all these ratios, we adopt an 
experimental point of view. We will choose a particular input boundary condition 
scenario and particular values of mi/2 and tan/? as 'nature's choice'. The resulting 
model will predict certain m^o and m~± values, which will be measured with small 
errors. The same values for these two observable masses can only be obtained for 
very specific mi/2 and tan (3 values in any other boundary condition scenario. Once, 
the (mi/2,tan/3) location in each scenario that yields the observed m^o and m~± 
is established, we compute the predictions for all the ratios of branching fractions. 
We use the notation Ri, with i specifying any particular ratio; the values of the 
Ri for the input scenario will be denoted by R^. We also compute the la error 
in the measurement of each of these ratios (denoted ARi) as found assuming that 
the input model is nature's choice. We may then compute the expected Ax^ for 
any of the other models relative to the input model as: 

Ax' = E Ax? , with Ax^ = . (11) 

We will see that very large Ax^ values are typically associated with an incorrect 
choice of model. 

It is important to note that many other observables that discriminate between 
models will be available from other experimental observations. An additional Ax^ 
contribution should be added for each observable in assessing the overall improb- 
ability of a model other than the correct one. However, there are advantages to 
restricting oneself to the branching fraction ratios only. For example, (which 
will be readily measured in slepton pair production) differs substantially at fixed 
m^o,mi/2 as one moves between the NS, D and HS scenarios, and would readily 
distinguish between the models. However, mj^ is primarily sensitive to the value 
of the slepton mo at M^/, which could differ from the mo associated with the Higgs 
fields if the GUT boundary conditions are nonuniversal. In contrast, the branch- 
ing fraction ratios are primarily sensitive to the Higgs mo value relative to mi/2. 
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Different sets of observables will have maximial sensitivity to different subsets of 
the GUT scale boundary conditions. The Higgs branching fraction ratios should 
be very valuable in sorting out the correct relation between mi/2 and the mo for 
the Higgs fields, and in determining tan /3. 

4.4 A Test Case 

As a specific example, suppose the correct model is D~ with mi/2 = 201.7 GeV 
and tan 3 = 7.50. This would imply m^o = 349.7 GeV, m~± = 149.5 GeV. The 
mi/2 and tan/5 values required in order to reproduce these same m^o and m~± 
values in the other scenarios are listed in Table |l[ Also given in this table are 
the predicted values of uiho and mj^ for each scenario. In order to get a first 
feeling for event numbers and for the errors that might be expected for the ratios 
of interest, we give in Table ^the numbers of events. A/" and T>, predicted in each 
scenario for use in determining the numerators and denominators of Eqs. (0)-® 
and Eqs. (§)-(§), assuming Leg = 80 fb~^ at ^/s = 1 TeV. These numbers include 
the SUSY branching fractions, B^s of Eq. (||), and so forth following the itemized 
list of factors given earlier 

Table 1: We tabulate the values of mi/2 (in GeV) and tan /3 required in each of our 
six scenarios in order that m^o = 349.7 GeV and m~± = 149.5 GeV. Also given 
are the corresponding values of 711^0 and mj . Masses are in GeV. 







D+ 


NS- 


NS+ 


HS- 


HS+ 


mi/2 
tan /? 
niHO 
mj 


201.7 
7.50 
350.3 
146.7 


174.4 
2.94 
355.8 
127.5 


210.6 
3.24 
353.9 
91.0 


168.2 
2.04 
359.0 
73.9 


203.9 
12.06 
350.1 
222.9 


180.0 
3.83 
353.2 
197.4 



From Table |^, we observe that the event rates for the /i > scenarios 

are all rather small as compared to the event rates for the /i < scenarios. (This 
happens because the mi/2 and tan/3 values required for m^o = 349.7 GeV and 
m~± = 149.5 GeV when yU > are very close to the scenario boundary.) For 
example, if the D~ model is nature's choice, the if°y4°-pair denominator rates 
would be ~ 198, implying a statistical error of only ~ ±14. Assuming systematic 
error of order 10%, the net error in event number would certainly be < 35, i.e. many 
a away from any of the /i > scenario predictions. We also see significantly larger 

''Because B{A^ tt) ~ and B{H^ tt) is typically small for the test case choice of 
TO^o = 349.7 GeV (given nit — 175 GeV), the ratio of Eq. (||) and its numerator event rate are 
both smaU. Note that Bcff{A° ^bb + tt) = B{A° bb) and that Boff ~>bb + tt) is not very 
different from B{H'^ bb) for this same reason. 
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Table 2: We give the numbers of events predicted in each scenario at the parameter 
space locations specified in Table |I] available for determining the numerators and 
denominators of Eqs. (|)-(|) and Eqs. These event rates are those for 

Leff = 80 fb~^ at \/s = 1 TeV. They include all branching fractions. Our notation 
is and ^^(#) for the event rates in the numerator and denominator, respectively, 
of the ratio defined in Eq. (#). 





D- 


D+ 


NS- 


NS+ 


HS- 


HS+ 




97.0 


92.3 


88.3 


49.2 


76.1 


124.0 


^§ 


0.1 


0.7 


3.8 


1.02 


0.0 


0.2 




16.4 


2.7 


46.6 


1.47 


3.8 


2.4 




2.0 


1.3 


9.2 


0.6 


0.4 


1.1 




198 


9.6 


62.1 


2.6 


250 


18.2 




198 


8.9 


58.3 


1.6 


250 


18.0 




225 


189 


138 


135 


189 


262 


^§ 


58.4 


4.2 


6.5 


1.1 


90.0 


9.5 


^§ 


13.0 


12.8 


21.9 


9.0 


3.3 


12.3 




317 


415 


445 


465 


320 


348 



numerator rates N"^ and A/(|) for the /i < scenarios than for the /i > scenarios. 
Thus, in this particular case, even before examining the branching fraction ratios, 
the /i > scenarios could be excluded. 

The Af and T> event numbers of Table ^ also make apparent the accuracy with 
which the ratios of Eqs. (0)-® and Eqs. (|])-(|) can be measured. For example, 
the event numbers A/| and show that good statistical precision, ~ ±10% — 15%, 
can be expected for the ratio of Eq. in the /i < scenarios. Such statistical 
precision implies that this ratio will also clearly distinguish between the input D~ 
scenario and any of the /i > model predictions. 

To illustrate the value of the branching fraction ratios more clearly, we present 
in Fig. |T9| a plot which gives the expected values and the ilex errors as a function 
of scenario for four of the ratios that will be useful in distinguishing between the 
different scenarios at the given input (measured) values of m^o and ttl~±. In this 
plot, the errors as a function of scenario are those that are expected if the scenario 
listed on the horizontal axis is the correct one. Thus, if the correct model is D~, 
the central value and ±la upper and lower limits for each ratio are those given 
above the D~ scenario label on the x-axis. The ability of each ratio to discriminate 
between a given scenario on the horizontal axis and one of the five alternatives is 
indicated by the extent to which the ±la error bars for the given scenario do not 
overlap the central points for the other scenario. Referring to Fig. ^ we observe 
the following. 
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Scenario Overlap of Branching Fraction Ratios 



E_ = 1 TeV, Leff = 80 fb"' 
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Figure 19: We plot the branching fraction ratios B{H'^ — > t'^v) / B{H~^ — 
th), B{H+ SUSY [Oi][> Oj])/B{H+ tb), B{H^,A^ 
SUSY)/B^s{H^\A^ bb,tt) and B{H+ h^W+)/B{H+ tb) with 
±1(7 error bars as a function of scenario, adjusting aiid tan/3 in each 
scenario so that m^o = 349.7 GeV and m~± = 149.5 GeV are held fixed. 

Error bars are for Leff = 80 fb~^ at i/s = 1 TeV, and are those that would 
arise if the input (nature's choice) scenario is that listed on the horizontal 
axis. No error bar is shown for the T^v/th ratio in the NS^ scenario since 
the predicted rate is less than 4 events; a very large error bar should be 
assumed. 
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• The ratio B{H+ SUSY [Oi][> Oj])/B{H+ th) succeeds in distin- 
guishing the D~ scenario from all but the HS^ and HS"*" scenarios. 

• The ratio B{H~^ t^v) / B{H^ tbh) provides excellent discrimination 
between the D~ input scenario and the D"*", NS~, NS"*", and HS"*" scenarios, all 
of which must have tan /3 < 4 (in order to reproduce m^o = 349.7 GeV, m~± = 
149.5 GeV) as compared to tan 13 = 7.5 for the D~ scenario. The much smaller 
tan/? values imply much smaller r^v/ih ratios, as was illustrated in Fig. |16|. 
The more limited ability of this ratio to discriminate between the high tan (3 
values of 7.5 for D~ vs. 12 for HS^ is also apparent from Fig. 0. 

• The ratio B{H°,A^ SU SY) / B{H^ , bb,tt) will strongly rule out 
fi > scenarios if ;U < is nature's choice. Due to the small error bars, this 
ratio provides some discrimination between the D~ and HS~ scenarios even 
though the predicted central values are not very different. 

• The ratio B{H+ h^W+)/B{H+ th) is quite different for the D", NS", 
and HS"*" scenarios as compared to the D+, NS"*", and HS~ scenarios. However, 
discrimination power is limited by the relatively large error bars. Nonetheless, 
this ratio yields a bit more than 2.1a discrimination against the HS~ model 
if the D~ model is nature's choice. 

The quite substantial dependence of the ratios on scenario and location in pa- 
rameter space, as displayed in Figs. |^, |1^, |ll|, |12|, |13|, |14|, and |1^, suggests that 



similar discrimination will be possible for most input scenario and parameter space 
location choices. 

In Table |^ we more thoroughly quantify the process of excluding the D+, NS~, 
HS~, and HS"*" scenarios relative to the input D~ scenario. There we give the 
contribution to Ax^ (computed relative to the assumed-to-be-correct D~ scenario) 
for each of a selection of independently measurable ratios. Also given for each of 
the incorrect scenarios is the sum of these contributions. This table shows that the 
D~ scenario can be distinguished from the D^, NS~, NS^, and HS^ scenarios at 
an extremely high statistical level. Further, even though no one of the branching 
fraction ratios provides an absolutely clear discrimination between the D~ and the 
HS~ scenarios, the accumulated discrimination power obtained by considering all 
the ratios is very substantial. In particular, although the ratios of Eq. (^), (H), and 
(^) are only poorly measured for Leg = 80 fb~^, their accumulated Ax^ weight can 
be an important component in determining the likelihood of a given model and 
thereby ruling out incorrect model choices. 

Thus, consistency of all the ratios with one another and with the measured 
m^o, neutralino and chargino masses will generally restrict the allowed models to 
ones that are very closely related. The likelihood or probability associated with 
the best fit to all these observables in a model that differs significantly from the 
correct model would be very small. 
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Table 3: We tabulate Ax?, see Eq. (|TT|), (relat ive to the D scenario) for the in- 
dicated branching fraction ratios as a function of scenario, assuming the measured 
m^o and m~± values are 349.7 GeV and 149.5 GeV, respectively. The SUSY chan- 
nels have been resolved into final states involving a fixed number of leptons. The 
error used in calculating each Ax? is the approximate la error (as defined in text) 
with which the given ratio Ri could be measured for L^g = 80 fb~^ at y/s = 1 TeV 
assuming that the D~ scenario is the correct one. 



Ratio 


D- 


D+ 




NS+ 


HS- 


HS+ 


2{B{H'\A" ^ SUSY ^ [Ol][> Oj])/ 





12878 


1277 


25243 


0.77 


10331 


B,s{H",A° ^bb,tt)) 














2{B{H^,A^ SUSY [U][> Oj])/ 





13081 


2.41 


5130 


3.6 


4783 


Bcs{H°,A° bb,tt)) 














2{B{H^,A^ SUSY [2i][> Oj])/ 





4543 


5.12 


92395 


26.6 


116 


B,s{H°,A° ^bb,tt)) 














B{H^ h°h°)/B{H° bh) 





109 


1130 


1516 


10.2 


6.2 


2B{H+ SUSY [0£][> Oj])/ 





12.2 


36.5 


43.2 


0.04 


0.2 


B{H+ tb) 














2B{H+ SUSY [!£][> Oj])/ 





1.5 


0.3 


0.1 


5.6 


0.06 


B{H+ tb) 














2B{H+ h"W)/B{H+ ^ tb) 





0.8 


0.5 


3.6 


7.3 


0.3 


2B{H+ Tv)/B{H+ tb) 





43.7 


41.5 


47.7 


13.7 


35.5 







30669 


2493 


124379 


68 


15272 



4.5 Separating Different SUSY Decay Modes 

An important issue is the extent to which one can be sensitive to the branching 
fractions for different types of SUSY decays of the Higgs bosons, relative to one 
another and relative to the overall SUSY decay branching fraction. Interesting 
SUSY decay rates include: 

• B{H^, —>■ XiXi + '^'^)) leading to a totally invisible final state; 

• B{H^, AO where i^Xi or uxf; 

• B{H^, A^ xtxl), where xf ^ ^^^xl, 33x\ or i^u (with i^V ^ £±xVx?); 

• B{H^ I^V), where i^Xi, or uxf; 

• B{H^ -> xtxi), where xf i^uXi, JJX? or i^u (with I^u i^Xi^^Xi)- 

Predictions for such rates depend in a rather detailed fashion upon the SUSY 
parameters and would provide valuable information regarding the SUSY scenario. 
For example, in going from NS to D to HS the masses of the sneutrinos and sleptons 
increase relative to those for the charginos and neutralinos. The ^4° 
and i^V branching fractions should decline in comparison to H^, A^ 
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XiXi ^-iid H"^ xfxi^ respectively. In small sections of the D and NS scenario 
parameter spaces, the sleptons and sneutrinos are sufficiently light that Xi decays 
almost exclusively to £^z/ followed by — »• i'^Xii^Xi^ implying that Xi decays 
would mainly yield leptons and not jets. 

The difficulty is that several different SUSY channels can contribute to any 
given final state. Two examples were noted earlier: the + I^t channel receives 
contributions from both A° and xtXi decays; and the i'^+I^T channel 

receives contributions from i^u and XiXi- Another example, is the purely 

invisible or final state; it can arise from either XiXi o^' ^"^^ (with u — >■ 
I'Xi) production. Thus, the physically distinct channels, defined by the number 
of leptons and jets present,^ typically have multiple sources. Still, a comparison 
between the rates for the final states so-defined might be quite revealing. For 
instance, if xf ~^ is not kinematically allowed, the XiXi final states are 
expected to yield more li + 2j and 0£ + 4j events than 2i + 0j events, whereas 
events will yield only 2i+0j events. Further, the ts must be of the same type in this 
latter case. The effective branching fraction for XiXi ~^ i'^i'+J^ t with both £'s of 
the same type is only 1/81. In addition, the £'s in the latter derive from three-body 
decays of the Xii ^^nd would be much softer on average than £'s from Even 
if this difference is difficult to see directly via distributions, it will lead to higher 
efficiency for picking up the events. Of course, if event numbers are sufficiently 
large (which in general they are not) that detailed kinematical distributions within 
each final state could be obtained, they would provide additional information. We 
do not pursue this latter possibility here. 

Based on the above discussion, the following ratios would appear to be poten- 
tially useful. For the and we consider: 



5(/70 


hh)B{A'' ■ 


+ B{A^ - 


bb)B{H^ - 




(12) 




hh)B{A^ 


SUSY) +5(^0 - 


bb)B{H^ - 


-^SUSY) ' 


B{H^ 


bb)B{A^ - 


^ [2i][0j]) + B{A' - 


bb)B{H^ - 


-mm) . 


(13) 




^ bb)B{AO 


SUSY) +5(^0 - 


bb)B{H^ - 


-^SUSY) ' 


B{H° 


bb)B{A^ 


[> Oim) + B{A'- 


bb)B{H^ - 




(14) 


B{HO 


bb)B{A^ 


SUSY) +5(^0 - 


bb)B{H^ - 


SUSY) 


BiH"^ ^ 


bb)B{A^ 


[Oi][> lj])+B{A'- 


bb)B{H^ - 


.m[> ij]) 


(15) 


B{HO 


bb)B{AO 


SUSY) +5(^0 - 


bb)B{HO - 


SUSY) 


5(i7° 


bb)B{A^ 


[li][> lj])+B{A'- 


bb)B{H^ - 


-m> ij]) . 


(16) 


B{HO 


bb)B{A^ 


SUSY) +5(^0 - 


bb)B{H^ - 


SUSY) 



(As before, rriAO rriHO implies that we cannot separate the and A^ via the 
tagging procedure.0) Once again, we employ shorthand notations for the quantities 

^The totally invisible final state would be [0-^][0j], and so forth. 

^The A" — > £^£^ branching ratio turns out to be rather small in the three GUT scenarios 
studied — the required L — R mixing is numerically very small in the slepton sector. 
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appearing in Eqs. (p!^)-([T^). For example, the ratio of Eq. (|T^) will be denoted by 

(17) 



{B{A'^,H^ ^ [0£][Oj])E(ij'°, A" -> bb)) B{H\A^ [Of][Oj]) 

{B{A\ SUSY)5(ifo, AO ^ 65)) 



B{H^,A^ SUSY) 



eff 



in what follows. 

For the we consider the ratios: 



B{H+ 




^6t)H 


h B{H- - 




^tb) _ 


B{H+ 


S\JSY)B{H- 


^6t)H 


h B{H- - 


SVSY)B{H+ ^ 


tb) ' 


B{H+ 


[> Um)BiH- 


^6t)H 


h B{H- - 


Um)B{H^ 


- ^ t6) 


B{H+ 


SVSY)B{H- 


^6t)H 


h B{H- - 


SVSY)B{H+ ^ 


tb) 


B{H+ 


[Oi][> lj])B{H- 


^bt)- 


VB{H- - 


^m[> l3])B{H^ 


^ ^ t6) 


B{m 


SVSY)B{H- 


bt) + B{H- - 


SVSY)B{H+ - 


>t6) 



(19) 



(20) 



The ratios of Eqs. (|18D-(|0D reduce to: 

B{H+ [l£][Oj]) B{H+ ^ [> U]m B{H^ 



> IJ]) 



5(iJ+ ^ SUSY) ' 5(iJ+ ^ SUSY) 



B{H+ SUSY) 



(21) 



respectively. 

Also of interest are ratios of the different numerator terms to one another 
within the above neutral and charged Higgs boson sets. All the ratios that one 
can form have the potential to provide important tests of the Higgs decays to the 
supersymmetric particle pair final states. 



To illustrate, we present two figures. In Fig. |2y we present three-dimensional 
lego plots of the ratio of Eq. ( |T^ ) as a function of location in (mi/2,tan/5) param- 
eter space. (Because of the combination of slow variation and very sharp changes, 
the contour plots similar to those presented earlier are rather difficult to inter- 
pret.) In Fig. we plot the numerator of Eq. (|^) divided by the numerator 
of Eq. (|20|) . In both sets of lego plots, the ratio is set to zero if there are fewer 
than 4 events in the numerator or denominator after including the earlier-discussed 
tagging/reconstruction efficiencies and assuming y/s = 1 TeV and Lcs = 80 fb~^. 

The most important feature apparent from these figures is the generally de- 
creasing magnitude of these two ratios as one moves from the NS to the D to the 
HS scenario. This is a result of the decreasing importance of slepton/sneutrino- 
related decays as compared to chargino/neutralino-based decays. When the lat- 
ter types of decay are prevalent, a much larger fraction of the events will have 
jets than if the former decays dominate. The decreasing importance of the slep- 
ton/sneutrino class is to be expected due to the increasing mass of these states as 
mo increases in going from NS to D to HS. The occasionally very large values of 
B{H+ [> li][Oj])/B{H+ [Oi][> Ij]) in Fig. |l|injhe D" and NS+ plots occur 
in the small wedges of parameter space where Xi ^^i^ decays are kinematically 
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Figure 20: We present lego plots of the ratio of Eq. (|T4|) in each of the six 
scenarios as a function of location in (772^/2 > tan /3) parameter space. The 
value of the ratio is given by the height on the z-axis. Non-zero values of 
the ratio are given only in regions where there are at least 4 events in the 
numerator after including tagging/reconstruction efficiencies. 
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B(H^ ^ SUSY ^ m 1 l][Oj])/B(H^ ^ SUSY ^ [0I][^ 1 jl) 

eV^H-'H', E,„=1 TeV, U„=80 fb'' 




Figure 21: We present lego plots of the numerator of Eq. ( p^ ) divided 
by the numerator of Eq. (|20| ) in each of the six scenarios as a function of 
location in (mi/2;tan/3) parameter space. The value of the ratio is given 
by the height on the z-axis. Non-zero values of the ratio are given only 
in regions where there are at least 4 events in both the numerator and 
denominator after including tagging/reconstruction efficiencies. 
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allowed, and final states containing only jets must arise from higher ino states and, 
thus, are very rare. 

It should be apparent from these two figures that rather dramatic differences 
between the scenarios at a given (mi/2,tan/5) location are the norm. In general, 
statistics are such that the different scenarios can be distinguished from one another 
at a substantial level of significance just on the basis of these two ratios. Ratios 
other than the two plotted ones can also provide good discrimination. We shall 
illustrate this for our standard m^o = 349.7 GeV, m~± = 149.5 GeV point discussed 
in association with Tables |I| and 0, Fig. jTU] and Table |^. 

Table m gives the {mi/2, tan j3) parameters required for m^o = 349.7 GeV, m~± = 
149.5 GeV in each of the six GUT scenarios. In Table § the event rates for the 
SUSY final states corresponding to the numerators of the ratios listed in Eqs. (|12|)- 
([TBD and (|lD-(|3) given for these (mi/2,tan/?) values. We will follow the same 
notation in terms of A/(Eq. #) as for Table |^. An examination of Table ^ reveals 
event rates in the individual channels that vary from a few events, implying poor 
statistics, to 50 or 60 events, for which statistical accuracy would be quite reason- 
able. 



Table 4: For the (mi/2,tan/3) values required for uiao = 349.7 GeV, m~± = 
149.5 GeV, we tabulate the numbers of events predicted in each scenario in the 
final states corresponding to the numerators and denominators of Eqs. (|T^-(rB) 
and ([r8| ) - (|20|) . These rates are those obtained for L^s = 80 fb~^ at y/s = 1 TeV. 
They include all branching fractions. 
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Not surprisingly, the ratios of rates of the various SUSY channels can contribute 
significantly to our ability to discriminate between different GUT scenarios. To 
illustrate, we follow the same procedure as in Table ^. Taking m^o = 349.7 GeV 
and m~± = 149.5 GeV, we assume that the correct scenario is D~ and compute 
the Ax^ by which the prediction for a given ratio in the other scenarios deviates 
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from the D~ prediction. Statistics are computed on the basis of the expected D~ 
rates, as given in Table |^. The resuhing Ax^ values are given in Table ^. Since 
these ratios are not all statistically independent of one another, we do not sum 
their Ax^'s to obtain an overall discrimination level. However, a rough indication 
of the level at which any given scenario can be ruled out relative to the D~ is 
obtained if we add the largest Axf from the neutral Higgs list and the largest from 
the charged Higgs list. The weakest discrimination level following this procedure 
is Ax^ ~ 15 in the case of the D"*" scenario. Note that this scenario is highly 
unlikely on the basis of the earlier J2i Axj value listed in Table 0. In Table ^, 
the weakest discrimination was that for the HS~ scenario with J2i ^Xi ~ 68. We 
observe from Table | that the ratio B{H^,A^ [0£] [0j])/S(i7°, A° -> [2£][0j])|cff 
has Ax^ ~ 928 for the HS^ case, which would certainly rule it out. 

Table 5: We tabulate Ax^, see Eq. (|llD, (relative to the D~ scenario) for the 
indicated ratios as a function of scenario, assuming the measured m^o and m~± 
values are 349.7 GeV and 149.5 GeV, respectively. The SUSY channels have been 
resolved into final states involving a restricted number of leptons and jets. Only 
those ratios with substantial power for discriminating between scenarios are tab- 
ulated. The error used in calculating each Ax^ is the approximate la error (as 
defined in text) with which the given ratio Ri could be measured for Leg = 80 fb^^ 
at ^/s = 1 TeV assuming that the D~ scenario is the correct one. 
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The above illustrations demonstrate that the ratios of rates for individual SUSY 
channels correlate strongly with the underlying physics of the different GUT sce- 
narios (light vs. heavy sleptons in particular) and add a powerful component to 
our ability to determine the correct scenario. 

5 Summary and Conclusions 

In this paper, we have considered detecting and studying the heavy Higgs 
bosons of the minimal supersymmetric model when pair produced in e"'"e~ or 
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collisions. We have shown that, in the SUSY GUT models studied, the target lumi- 
nosities of L = 200 fb~^ and L = 1000 fb~^ at y/s = 1 TeV and 4 TeV, respectively, 
will allow detection of H^A^ and H^H^ pair production throughout essentially all 
of the model parameter space which is allowed by theoretical and kinematic con- 
straints, despite the presence of SUSY decay modes of the , A^, at a signifi- 
cant level. The all-jet and high-multiplicity final states coming from H^, A^ — > bb, tt 
and — >■ tb, H~ — >■ 6t are essentially background free and provide appropriate 
and efficient signals with rates that are adequate even when SUSY decays are 
present. In the all-jet channels, the individual Higgs boson masses, m^o, muo and 
m//+, can be measured and the approximate degeneracy (m^o ~ 111^0 ~ rnH±) 
predicted by the MSSM can be checked. 

Once the Higgs bosons are detected and their masses determined, the relative 
branching fractions for the decay of a single Higgs boson can be measured by 'tag- 
ging' [i.e. identifying) one member of the H^A^ or H^H'^ pair in an all-jet mode, 
and then looking at the ratios of the numbers of events in different event classes 
on the opposing side. In this way, the relative branching ratios of Eqs. 
Eqs. (||)-(||), Eqs. (|T^)-(p!6D, and Eqs. ([T8|)-(p!5D can be measured with reasonable 
accuracy whenever parameters are such that the final states in the numerator and 
denominator both have significant event rate.j] We find that the measured Higgs 
masses and relative branching fractions, in combination with direct measurements 
of the chargino and neutralino masses, will over-constrain and very strongly limit 
the possible SUSY GUT models. 

The specific SUSY GUT models considered are moderately conservative in that 
they are characterized by universal boundary conditions. In all, we delineated ex- 
pectations for six different models, requiring correct electroweak symmetry break- 
ing via evolution from the GUT scale to uiz- For each model, there are only two 
parameters: mi/2 (the universal gaugino) mass; and tan/3 (the usual Higgs field 
vacuum expectation value ratio). Each model is characterized by a definite relation 
of the universal soft-SUSY-breaking scalar mass, mo, and the universal mixing pa- 
rameter, Ao, to mi/2, as well as by a choice for the sign of /i (the Higgs superfield 
mixing coefficient). 

The strategy for checking the consistency of a given GUT hypothesis is straight- 
forward. First, the measured ^4°, neutralino and chargino masses are, in almost 
all cases, already sufficient to determine the mi/2 and tan/3 values required in the 
given GUT scenario with good precision. The value of tan (3 so obtained should 
agree with that determined from chargino pair production rates. The Higgs sec- 
tor branching fractions can then be predicted and become an important testing 
ground for the consistency of the proposed GUT hypothesis as well as for testing 
the MSSM two-doublet Higgs sector structure per se. 



^We focus on event rate ratios rather than the absolute rates in the many different channels 
since the possibly large systematic errors of the absolute rates will tend to cancel in the ratios. 
In some cases, absolute event rates are so different that they would also provide substantial 
discrimination between different models, despite the possibly large systematic errors. 
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Within the hst of ratios of branching fractions given in Eqs. (@)-(^) and 
the averagel H^, SUSY, the H+ SUSY and the H+ t+v branching 

fractions typically fix a relatively precise location in (mi/2, tan (3) parameter space. 
These values can be compared to those required by the m^o and m~± mass mea- 
surements. Consistency within experimental errors is typically only possible for 
a small set of closely related models. In the sample situation detailed in Section 
4, where we assumed that one of the six GUT models was correct and computed 
statistical errors on that basis, only one of the remaining five models could possi- 
bly be confused with the input model after measuring the above three branching 
fractions relative to that for the final state used for tagging. By subdividing the 
SUSY signal into final states with a definite number of leptons and any number 
of jets, and considering as well the H*^ — h^h^, Zh^ and W^h^ 

branching fractions, we found it possible to distinguish between these two choices 
at a very substantial statistical level. Thus, a unique model among the six rather 
similar models is singled out by combining measurements from the Higgs sector 
with those from conventional SUSY pair production. In short, measurements de- 
riving from pair production of Higgs particles can have a great impact upon our 
ability to experimentally determine the correct SUSY GUT model. 

The above discussion has left aside the fact that for universal soft-scalar masses 
the measured value of the slepton mass would determine the relative magnitude of 
mo and mi/2- Of the two models mentioned just above, one has a large mo/mi/2 
value and the other a much smaller value. They could be easily distinguished on 
the basis of mj alone. However, if the soft-scalar slepton mass is not the same as 
the soft-scalar Higgs field masses at the GUT scale, the branching fraction ratios 
would give the best indication of the relative size of the soft-scalar Higgs mass as 
compared to mi/2- 

More information regarding the slepton/sneutrino mass scale and additional 
ability to discrminate between models are both realized by subdividing the SUSY 
decays of the Higgs bosons in a way that is sensitive to the relative branching 
fractions for slepton/sneutrino vs. chargino/neutralino decays. Slepton/sneutrino 
channels essentially only produce leptons in the final state, whereas the jet compo- 
nent is typically larger than the leptonic component for chargino/neutralino decays 
(other than the totally invisible mode). Thus, we are able to define individual 
SUSY channels, characterized by a certain number of leptons and/or jets, which 
display a strong correlation with the slepton/sneutrino decay component. We find 
that these individual channels have sufficiently large event rates that the ratios of 
the branching fractions for these channels can typically be determined with reason- 
able statistical precision. For the earlier-mentioned input model, we can compute 
the statistical level at which the other five GUT scenarios would be ruled out us- 
ing various of these ratios of branching fractions. Excellent discrimination between 
models on this basis is found. 

^Only the average can be determined given that typicaUy m^o ^ m/^o. 



52 



In conclusion, our study shows that not only will detection of Higgs pair pro- 
duction in e+e~ or /i+zi" collisions (at planned luminosities) be possible for most of 
the kinematically accessible portion of parameter space in a typical GUT model, 
but also the detailed rates for and ratios of different neutral and charged Higgs 
decay final states will very strongly constrain the choice of GUT-scale boundary 
conditions. In estimating experimental sensitivity for Higgs pair detection and for 
measuring Higgs masses and branching fractions, we included substantial inefficien- 
cies and all relevant branching fractions. Although we believe that our estimates 
are relatively conservative, it will be important to re- visit this analysis using a full 
Monte Carlo detector simulation. 
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